Encapsulation and polymerization of acetylene molecules inside a carbon nanotube 
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We study the energetics of acetylene (C2H2) molecules inside a carbon nanotube (CNT) using 
the ab initio pseudopotential method. The encapsulation energy of a single C2H2 molecule into 
the nanotube and the formation energy of (C2H2) n @CNT are calculated. We investigate whether 
a polyacetylene chain can be produced by fusion of the close-packed acetylene molecules inside the 
CNT and find that there is practically no activation barrier to polymerization. We propose to 
employ this method to obtain straight, perfectly isolated, and single-stranded polyacetylene chains 
encapsulated inside CNTs, which may be used for molecular electronic devices. 



PACS numbers: 61.46.+W, 73.22.-f, 82.35.Cd 

I. INTRODUCTION 

Due to the presence of the inner hollow space, the CNT 
can be an ideal container for atoms or small molecules. 
Recently, Takenobu and coworkers reported that organic 
molecules predominantly occupy the inner space of the 
nanotubesi Atoms and molecules encaged in the nan- 
otube are isolated from the environment and their chem- 
ical reaction could be triggered by electron beam irradi- 
ation. Among alkynes (organic molecules with carbon 
triple bonds), acetylene (C2H2) is the simplest and a 
good model system to study pressure-induced polymer- 
ization which arises at pressures above 3 GPa at room 
temperature^. On the other hand, under pressure or 
in the presence of copper, it can be converted to car- 
bon and hydrogen with explosion. This instability causes 
problems in the handling and storage of the material and 
the commercial acetylene is stabilized with acetone and 
phosphine to prevent explosion. 

Mcintosh and coworkers proposed to insert the poly- 
acetylene (PA) chains inside carbon nanotubes and avoid 
interchain coupling and complicated morphology of the 
bulk materials^. Because of the entropic effect, how- 
ever, the PA (polymer) is not straight and cannot reach 
the deep inside of the CNT at room temperature. Here 
we propose a method to store acetylene molecules and 
convert them into PA using the CNT as a chemical reac- 
tor as well as a container. Through this process, we can 
obtain a straight, perfectly isolated, and single-stranded 
PA encapsulated inside a CNT. Carbon fullerenes were 
already found to be converted into a CNT inside the host 
nanotube by the coalescence process§i A collection of 
PA chains under such a morphological control may be 
used for future molecular electronic devices. In this work, 
we calculate the energetics of acetylene insertion into a 
nanotube and find that the reaction is exothermic result- 
ing in a static capillary forced. This force compresses 
incorporated acetylene molecules with an effective pres- 
sure of the order of sub-GPa. Once acetylene molecules 
are aligned compactly inside the tube, polymerization of 
them is predicted to take place practically without an 
activation energy barrier. 



II. COMPUTATIONAL DETAILS 

We perform ab initio electronic structure calcula- 
tions for C2H2 molecules inside the CNT based on 
the density functional theory (DFT) within the Local 
Density Approximation for the exchange-correlation en- 
ergy. We also checked the energy shift by employ- 
ing the Generalized Gradient Approximation for some 
important geometries, and found that overall results 
are the same. For the large supercell calculations, 
norm-conserving pseudopotentials^ are employed in the 
Kleinman-Bylander formAfi. The wave functions are ex- 
panded using the numerical atomic orbital basis selpiiA^ 
and an energy cutoff of 100 Ry. To calculate the acti- 
vation energy barrier in polymerization of the acetylene 
molecules which require heavier computational loads, the 
nudged elastic band (NEB) methodi^i^ implemented in 
the PWSCF package^ is used. Ionic cores are described 
by ultrasoft pseudopotentials and the Kohn-Sham or- 
bitals are expanded in plane waves with a kinetic energy 
cutoff of 30 Ry. We have done test calculations for dia- 
mond and Ceo, and obtained converged results for elec- 
tronic and structural properties with two cutoff energies, 
100 and 30 Ry, for norm-conserving and ultrasoft pseu- 
dopotentials, respectively. The NEB method enables us 
to find the optimal reaction pathway with the minimum 
energy when both the initial and final states are known. 
Seven replicas are chosen including the initial and final 
ones to construct the elastic band in this method. The 
atomic positions of the CNT, acetylene, and polyacety- 
lene chain are relaxed, respectively until the atomic force 
becomes less than 0.05 eV/A. 



III. RESULTS AND DISCUSSION 

We first discuss the energetics of C2H2 insertion into 
a CNT. The dangling bonds of the open (5,5) nanotube 
of 2 nm in length are passivated by hydrogen atoms and 
the C2H2 molecular axis is assumed to be initially aligned 
with the tube axis as shown in Fig. 1(a). The supercell 
size in the lateral (xy) direction is chosen to be as large 
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as 20 A to avoid interaction between (5,5) CNTs with 
~ 7 A in diameter and that in the axial direction (z- 
axis) is 70 A to allow for sufficient vacuum regions. The 
abscissa in Fig. 1(b) represents the position of the center 
of C2H2 from the edge carbon atoms of the CNT along 
the tube axis. A static capillary force associated with 
the energy gain during encapsulation may compress other 
C2H2 molecules inside the CNT. The energy lowering by 
the axial insertion of an acetylene molecule into the (5,5) 
CNT is ps 0.65 eV as shown in Fig. 1(b), indicating an 
exothermic reaction. The energy slope along the tube 
axis in Fig. 1(b) corresponds to a static capillary force of 
ps 0.22 nN. The energy is further lowered by 0.1 eV when 
C2H2 is tilted by about 45° with respect to the axis of the 
(5,5) nanotube. In this situation, the distance between a 
H atom on C2H2 and the tube wall is ~ 2.3 A. Because 
of the limitation of the computational resource, the data 
presented in Fig. 1(b) are obtained without optimization 
for the whole system. According to our test calculation, 
the full relaxation lowers the total energy by only 0.03 
eV (30 meV). Thus we believe that we can do without 
complete optimization for the supercell containing the 
whole system (i.e., the nanotube and acetylene molecule). 

Small but finite energy variation for positive values 
of the abscissa in Fig. 1(b) implies that the acety- 
lene molecule inside the tube feels the lattice struc- 
ture of the CNT. The periodicity of the variation is 
almost the same as that of the armchair CNT (~ 
2.5 A). In addition, the tilting of the encapsulated 
C2H2 molecule indicates the nonbonding interaction be- 
tween the nanotube wall and C2H2. The recent den- 
sity functional calculations show that the long organic 
molecules such as tetracyanop-quinodimcthanc (TCNQ) 
and tetrafluorocyanop-quinodimcthane (F4TCNQ) en- 
capsulated in the bigger CNT are tilted by about 30° 
with respect to the tube axisi^iii. 

We then calculate the activation energy barrier along 
the reaction path between two locally stable configura- 
tions using the NEB method. Two stable configura- 
tions are the one-dimensional close-packed array of C2H2 
molecules (Geometry A) and the PA chain (Geometry B) 
inside the (5,5) CNT, respectively displayed in Fig. 2(a) 
and (b). As mentioned above, the acetylene molecule 
prefers to be tilted by 45° to the tube axis inside the (5,5) 
CNT (Geometry A). First of all, the formation energy of 
Geometry A is calculated. We define the formation en- 
ergy of the close-packed array of acetylene molecules in 
the (5,5) CNT, 

Ef orm = E[CNT + (C 2 H 2 )„] - E[CNT] - £[(C 2 H 2 )„], 

(1) 

where _B[CNT] and -E[(C2H 2 ) rl ] are the energies of the 
isolated nanotube and C2H2 while i?[CNT + (C 2 H 2 )„] is 
the energy of the whole encapsulated system. For the 
plane-wave calculation, the unit cell size in the lateral 
(xy) direction is chosen to be 15 A and that in the axial 
direction (z-axis) is 2.5 A with the periodic boundary 
condition. The formation energy is —0.4 eV per C2H2 
(exothermic reaction). We do not consider a Peierls dis- 
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FIG. 1: Energetics of acetylene insertion into a CNT. (a) 
Atomic model of a C2H2 molecule entering a (5,5) nanotube. 
(b) Energy profile during the insertion of C2H2 along the tube 
axis. 



tortion in calculating the energy of the PA chain. Figure 
2(c) reveals that Geometry B has a lower energy than 
Geometry A by 4.3 eV per C2H2 and there is no acti- 
vation barrier. In this reaction process, a 7r-electron in 
the triple bonding rapidly hops to form a new carbon- 
carbon bond with the nearest-neighbor C2H2. The initial 
and final dots in (c) indicate the energies of Geometries 
A and B, respectively. Therefore, close-packed acety- 
lene molecules are driven to be polymerized by fusion 
in the tube. During the process of polymerization, the 
carbon-carbon bonds of the acetylene molecule are elon- 
gated (from 1.2 to 1.4 A) and the angles between carbon 
and hydrogen atoms are changed. We have also tested 
the acetylene encapsulation and polymerization inside a 
larger diameter nanotube, (6,6) CNT. Insertion of C2H2 
into the (6,6) tube is more favorable than that of the 
(5,5) tube energetically. However, the (6,6) tube has a 
region of a constant potential inside^ and the random- 
ness caused by the lateral freedom of the C2H2 position 
in that region (owing to the entropic effect) makes poly- 
merization a little bit less favorable. For dense C2H2 
molecules in the (6,6) tube, the equilibrium angle of the 
C2H2 axis with the tube axis is found to be 90° before 
polymerization. Other than these differences, we obtain 
similar energetics to the (5,5) tube case, including the 
absence of the activation barrier to polymerization. 

Even if there is practically no activation energy barrier, 
it is very hard for the molecules to react with the neighbor 
molecule spontaneously inside the CNT before arrang- 
ing themselves into a one-dimensional close-packed array 
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FIG. 2: (color online) Polymerization process of acetylene in 
the CNT. (a) Atomic model of a one-dimensional close-packed 
array of acetylene molecules inside a (5,5) CNT (Geometry 
A), (b) Atomic model of the trans-polyacetylene formed inside 
a (5,5) CNT (Geometry B). (c) Energy diagram along the 
reaction pathway of polymerization from Geometry A to B. 

(Geometry A). The isomers of C4H4 cannot be created 
from acetylene molecules at room temperature because 
C2H2 molecules are not dissociated. In addition, some 
isomers of C4H4 such as cyclobutadien (ring-shaped) are 
energetically unstable. When the encapsulated C2H2 
molecules form the one-dimensional close-packed array 
by the effective pressure of the order of sub-GPa associ- 
ated with the capillary force, polymerization of acetylene 
takes place without activation energy barrier. Interest- 
ingly, it is known that the equilibrium spacing between 
Ceo molecules in the nanotube is smaller by ~ 4 % than 



in the three-dimensional Cqq crystal by the high resolu- 
tion transmission electron microscopy (HRTEM)A£, elec- 
tron diffraction 1 ^, and Raman measurements^. It indi- 
cates that incorporated molecules inside the CNT may be 
subject to an axial compressive force due to the energy 
change along the tube axis during the insertion^. In the 
same manner, acetylene molecules may be close-packed 
in the nanotube. 



IV. CONCLUSION 

We have performed ab initio pseudopotential calcula- 
tions to investigate the energetics of acetylene (C2H2) 
molecules inside a CNT. The encapsulation process of a 
single C2H2 molecule in the nanotube is an exothermic 
reaction and the static capillary force produced by enter- 
ing C2H2 is a fraction of a nanonewton. It is found that 
there is practically no activation barrier to the polymer- 
ization reaction, implying that a PA chain can be easily 
produced inside the CNT by fusion of the close-packed 
acetylene molecules in it. Thus, we propose a method to 
fabricate a straight and isolated PA chain inside a CNT, 
by which we can avoid the interchain coupling and under- 
stand the intrinsic properties of the individual PA chains. 
In addition, it is a promising way to investigate and con- 
trol the interaction between the CNT and the polymer. 
We believe that this study sheds light on physical and 
chemical properties of organic molecules and leads to an 
application of the nanotube as a nano-sized chemical re- 
actor. 
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